Human Mutation

CRISPR Cas9's ability to bind targeted and non-targeted sites and its functional interactions with other proteins and domains have been studied using many technologies. These studies are often limited by using oligos as substrates or by a reliance on predicted binding sites. We show how electronic genome mapping (EGM) with its ability to detect single-molecule DNA nicks can determine how gRNA mismatches impact targeted versus off-target binding across the whole human genome in vitro. Cas9 was found to be capable of binding a single site adjacent to the FXN repeat expansion responsible for Friedreich Ataxia using a gRNA specific for that sequence while gRNAs matching a repetitive sequence were able to bind hundreds of sites. Binding to perfectly matched gRNA sites was compared to sites with mismatches at each gRNA position. The PAM site and nearby seed region bases were confirmed as critical for binding. In contrast, single mismatches toward the gRNA 5' end allowed off-target binding. We also used EGM to determine how physically close the bound dCas9/gRNA complex can be relative to nearby DNA modification sites and still allow activity. Bound CRISPR/dCas9 inhibits nickases with recognition sites 4-7 bp away from the gRNA/PAM. Because EGM can simultaneously evaluate many mismatched off-target sites, these insights may be extended to other CRISPR systems and enable optimization of Cas proteins and gRNAs. This is the first example of EGM being used to characterize proteins and RNAs. This is also the first study demonstrating the combined application of CRISPR and EGM methodology and its potential role in the detection of specific repeat regions in the human genome that are associated with human repeat expansion disorders.

Our labs have been studying neurofibromin function and phenotype for over a decade with the intent of generating targeted therapeutics for Neurofibromatosis type 1 (NF1). In the process, we have generated numerous human cell lines containing variants within the NF1 gene. Herein, we present data characterizing these cell lines and make them publicly available for use by researchers both within and outside the NF1 community. We describe lines that contain both well-characterized patient-specific variants either at their endogenous locus or as exogenous cDNAs, as well as variants of uncertain significance (VUS), engineered as heterozygous, homozygous, and compound heterozygous variants. Methods to generate each line and subsequent validation steps are detailed including targeted sequencing, Western blot analysis for neurofibromin expression and ERK activation. The utility of each line is dependent on the variant of interest, the parental cell line, and the mechanism of action relevant to possible therapeutic targeting.

Background: Small nuclear RNAs (snRNAs) are RNA components of the major and minor spliceosomes that play a core role in splice-site recognition and control of the splicing process. Variants in genes that produce snRNAs are increasingly recognised as major contributors to rare disorders, including neurodevelopmental disorders (NDD) and retinal dystrophies (collectively termed RNUopathies, a subset of spliceosomopathies). Clinical interpretation of variants in snRNAs is, however, challenging and existing guidance to support clinical variant classification does not adequately capture the unique features of snRNAs that necessitate a bespoke approach. Methods: We quantified the elevated background mutation rate in snRNA genes using de novo variants from 12,007 trios and assessed mutation density in 76,215 genome sequenced individuals in gnomAD. We convened a panel of clinical, research, and industry scientists with wide-ranging expertise in clinical variant interpretation and classification and expert knowledge in snRNA genes to draft and refine a guidance document. Results: We detail important considerations for variant classification in snRNA genes. These include: the difficulties of variant identification which requires genome or targeted sequencing approaches, the large number of gene paralogs with high sequence identity that complicate read mapping and variant calling, and historical inaccuracies in snRNA gene annotation. Further we show a ~50-fold increase in de novo mutation rate in snRNA genes compared to intergenic sequence and discuss the implications of this for variant classification. We provide a set of specific recommendations for classifying variants in snRNA genes. Finally, we introduce RNUdb, an interactive web-based tool to support snRNA variant annotation and classification. Conclusions: We provide the first guidance for clinical variant classification in snRNA genes and anticipate that this will support routine screening and analysis of snRNA genes in clinical genetic testing.

Background: Short-rib polydactyly syndrome type IV (Beemer-Langer syndrome) is a lethal skeletal ciliopathy that can result from biallelic pathogenic variants in KIAA0586 (TALPID3). Earlier clinical observations suggested an unusually high occurrence of this phenotype in Romani families from the Carpathian Basin. Methods: Whole-exome sequencing was performed in affected individuals from five independent families, and the frequency of KIAA0586 c.1815G>A was assessed in a Hungarian Romani genetic database of 500 individuals. Results: The affected individuals were homozygous for c.1815G>A, a splice-disrupting variant previously shown to cause exon 14 skipping. Four families self-identified as Roma from Hungary and one affected individual was from Turkey. Six heterozygous carriers were identified among 500 individuals in the Hungarian Romani cohort, corresponding to an allele frequency of 0.6%. This was markedly higher than the frequencies cited from ALFA (0.01%) and gnomAD (0.0014%). The Hungarian-Romani and Turkish affected individuals also shared genotypes at informative markers flanking the variant, supporting a shared ancestral haplotype. Conclusions: These observations strengthen evidence that KIAA0586 c.1815G>A is a founder allele enriched in at least some Romani populations. The data support consideration of ancestry-informed carrier testing in appropriately counseled families. Historical admixture may explain the occurrence of the same haplotype in a Turkish case, but this remains a hypothesis requiring broader population-genomic confirmation.

FBXO31 encodes the substrate receptor of an SCF (SKP1-CUL1-F-box) E3 ubiquitin ligase that recognizes C-terminal amide-bearing proteins (CTAPs) and targets them for proteasomal degradation, a function of particular importance in neurons. A recurrent de novo missense variant, c.1000G>A p.(Asp334Asn), reprograms FBXO31 substrate specificity and causes an autosomal-dominant neurodevelopmental disorder with spastic cerebral palsy through a neomorphic gain-of-function mechanism. By contrast, biallelic loss of function has tentatively been linked to a recessive disorder. Here we report a female patient, born to consanguineous parents, with moderate intellectual disability, developmental delay and dysmorphic features, who carries a homozygous 54.3 kb deletion removing exons 2-8 and the coding part of exon 9 of FBXO31. The deletion retains only the exon 1-encoded F-box domain, while removing the coding sequence for the entire C-terminal substrate-recognition region. It is predicted to result in loss of function. Segregation is consistent with autosomal-recessive inheritance, with both parents and three of five unaffected siblings being heterozygous carriers. Together with two previously reported recessive families and the experimentally established neomorphic action of p.(Asp334Asn), the homozygous deletion supports an FBXO31 allelic series comprising dominant gain of function and recessive biallelic loss of function.

Lipase maturation factor 1 (LMF1) is a chaperone for lipoprotein lipase (LPL) and critically required for the enzyme to attain lipase activity. LMF1 has been identified as a canonical gene affected in Familial Chylomicronemia Syndrome (FCS) and it has also been extensively analyzed in Multifactorial Chylomicronemia Syndrome (MCS), the polygenic form of severe hypertriglyceridemia (hTG). While recent genetic studies resulted in an explosion in the number of hTG-associated LMF1 variants in different populations, their functional significance remains largely unknown. Here, we present a second-generation LMF1 activity assay allowing the streamlined quantitative functional analysis of LMF1 variants. The assay is based on the reconstitution of lipase maturation in transfected LMF1-deficient cells and fluorescence-based measurement of LPL activity in the culture medium. The use of Gaussia luciferase (GLuc)-LMF1 fusion constructs allows the assessment of LMF1 protein expression and the calculation of LMF1 specific activity. Several variants previously reported as neutral resulted in loss of LMF1 function in the new assay demonstrating increased sensitivity to detect variants of modest effect sizes. Furthermore, the functional analysis of 14 hTG-associated missense mutations identified 12 loss-of-function (LOF) variants. To gain initial insights into the structural determinants of LMF1 function, we analyzed the domain distribution of 73 hTG-associated missense variants and identified Loop C as a critical region in lipase maturation. In conclusion, we developed a simple, non-radioactive and sensitive quantitative assay of LMF1 activity to facilitate the functional analysis of hyperlipidemia-associated genetic variants and structure-function studies.

Non-coding variants altering mRNA splicing are increasingly recognized as important cause of Mendelian disorders. Deep intronic variants (DIV) that activate cryptic exons (CEs) cause ~20% inherited retinal diseases (IRD) cases, yet it remains challenging to predict intronic variant pathogenicity accurately, thus we used a high throughput splicing assay (HTSA) to measure the effects of rare deep intronic variants in recessive IRDs with one confirmed pathogenic variant. 640 very rare deep intronic variants were chosen in trans-position with known mutation from 76 patients. We tested them with HTSA, which consists of a split-GFP minigene, separated by an SMN1 gene intron into which the 270bp sequences flanking the DIVs were cloned. The plasmid library was transfected into 293HEK cells, and 48 hours later total RNA was extracted from the transfected cells. RNA transcripts produced by HTSA minigene were amplified by RT-PCR and sequenced. The intron sequences spliced between the two GFP exons were identified and quantified. 98 variants activated CE >100 times more in altered oligo compared to reference oligo and were classified as pathogenic in this experimental setting, validated in longer context, 26 variants were classified as VUS. Interestingly, only 6 from 90 variants (6,6%) were predicted to be pathogenic by in silico algorithms (Splice AI) which leaves room for significant improvement of prediction algorithms. The results of this experiment confirmed diagnosis for 50 of 78 patients (64%). So, 19,4% of deep intronic variants were shown to cause CE activation. This can lead to conclusion that significant portion of undiagnosed patients with recessive IRDs carry pathogenic intronic mutation which causes the disease.

CACNA1A-related disorders arise from altered function of the P/Q-type calcium channel Cav2.1 and encompass a broad spectrum of neurological phenotypes, including developmental and epileptic encephalopathy, familial hemiplegic migraine type 1 and episodic ataxia type 2. Although loss-of-function and gain-of-function mechanisms have both been implicated, interpretation of rare CACNA1A missense variants remains challenging because clinical presentation is heterogeneous and many variants identified through diagnostic sequencing remain classified as variants of uncertain significance. To generate clinically interpretable functional evidence, we developed a calibrated high-throughput automated patch-clamp assay to assess the CACNA1A variants. Twenty-six curated benign/likely benign and 30 curated pathogenic/likely pathogenic missense variants were expressed in Flp-In HEK293 cells with CACNB4 and CACNA2D1. We quantified current density, {Lambda}V0.5 activation,{Lambda}V0.5 inactivation and recovery from inactivation (T0.5 RFI), and used benign controls to define parameter-specific normal ranges with {+/-}3 Z-score thresholds. This multi-parameter framework classified all benign controls as functionally normal and 29 of 30 pathogenic controls as abnormal, yielding 100% specificity and 96.7% sensitivity. Current density identified 19 pathogenic controls with loss of function, while gating and recovery from inactivation identified additional abnormal variants with preserved or partially preserved current density. The assay generated likelihood ratios of 0.033 for normal function and 25.13 for abnormal function, supporting BS3 and PS3 functional evidence, respectively. Application to 23 CACNA1A variants identified in Australian patients showed abnormal function in all four variants already classified as pathogenic or likely pathogenic. Among 19 variants of uncertain significance, 11 were functionally abnormal and nine were reclassified as likely pathogenic after integration of PS3 evidence with other ACMG/AMP criteria. These data show that CACNA1A missense variants can produce loss-of-function, gain-of-function or mixed effects, and that multi-parameter functional assessment is required for accurate interpretation. This calibrated automated patch-clamp framework provides clinically interpretable evidence for refining CACNA1A variant classification in diagnostic practice.

BackgroundAcrodysostosis type 2 (ACRDYS2) is a rare autosomal dominant disease characterized by skeletal defects and cognitive deficit, with clinical symptoms observed in multiple other tissues including the skin. It is caused by mutations in a phosphodiesterase, PDE4D, a key regulator of cAMP/PKA (cyclic adenosine monophosphate / protein kinase A) signalling. Despite its well-defined genetic causes, the molecular mechanisms underlying the disease remain poorly understood, with studies based largely on engineered cellular models reaching conflicting interpretations. MethodsTo investigate how endogenous dynamics are affected by PDE4D mutations in unmanipulated cells, we studied PDE4D transcript and protein expression, activity and downstream signalling in native dermal fibroblast from ACRDYS2 patients and healthy controls. ResultsSignificant reduction in total PDE4D expression in patient cells was observed both at the transcript and protein level, with marked decreases in the long isoforms PDE4D4 and PDE4D7; a reduction in PDE4D9 mRNA was also observed. PDE4D enzymatic activity was reduced in ACRDYS2 fibroblasts, though total PDE activity was largely preserved. Reduced PDE4D expression was associated with an increase in the phosphorylated form of the cAMP-responsive transcription factor CREB and elevated PRKAR1A (PKA type 1 regulatory subunit alpha) transcript levels, suggesting altered downstream signalling. Interestingly, expression of the related phosphodiesterase family member PDE4B was increased, consistent with a compensatory response to reduced PDE4D function. ConclusionsThis is the first study demonstrating reduced PDE4D expression and isoform-specific dysregulation in native ACRDYS2 cells. Together, our results support a model in which reduction in PDE4D activity and compensatory changes in other PDE4 family members contribute to the molecular pathology of ACRDYS2, providing new insights into the molecular mechanisms underlying this disorder.

CBX1-related syndrome is characterized by developmental delay, hypotonia, autistic features, and mild dysmorphic features. This syndrome is caused by heterozygous missense variants in CBX1, which encodes heterochromatin protein 1 beta (HP1{beta}). These variants are situated within the chromodomain of HP1{beta}, a critical region that mediates the interaction between HP1{beta} and trimethylated H3K9. While a prior study identified three individuals with pathogenic CBX1 variants, the variant and clinical spectrum of CBX1-related syndrome remains largely unknown. Here we report three patients with novel de novo heterozygous variants in CBX1. Patients with CBX1 p.(Leu40Pro) and p.(Cys60Arg) manifested with global developmental delay, while the patient with CBX1 p.(Cys60Gly) primarily exhibited speech delay without concurrent motor developmental delay. All variants reside within chromodomain of the HP1{beta} protein. Interestingly, in vitro cellular assays revealed that CBX1 p.(Cys60Gly) variant showed an intermediate cellular phenotype between the wild type and the CBX1 p.(Leu40Pro) and p.(Cys60Arg) variants, suggesting a possible genotype-phenotype correlation. This report provides further evidence that missense variants within the chromodomain of HP1{beta} cause a syndromic neurodevelopmental disorder, CBX1-related syndrome.

Background: Short-read whole-genome sequencing (sWGS) has limitations in characterizing tandem-repeat expansions and variants in highly homologous or pseudogene-associated regions. Long-read sequencing can address some of these challenges but often benefits from high-molecular-weight DNA. TruPath Genome incorporates spatial proximity information into a short-read sequencing framework. This study evaluated its performance in characterizing selected complex genomic variants using archived clinical DNA specimens. Methods: Ten archived peripheral blood DNA specimens were analyzed, comprising seven cases with disease-associated short tandem repeat (STR) expansions and three cases involving homologous or pseudogene-associated loci (STRC, SMN1/SMN2, and CYP21A2). TruPath results were compared with previously available orthogonal assay results, including triplet repeat-primed PCR (TP-PCR), multiplex ligation-dependent probe amplification (MLPA), quantitative PCR (qPCR), and sWGS-based ExpansionHunter analysis. Results: TruPath generated repeat-length estimates across the evaluated STR loci. For two FMR1 full-mutation alleles, TruPath estimated 284 and 295 CGG repeats, respectively, while TP-PCR reported expansions exceeding 200 repeats. In a DMPK case, TruPath estimated an expanded allele of 1,203 CTG repeats, compared with a TP-PCR result exceeding 100 repeats. At the STRC and SMN1/SMN2 loci, TruPath findings were consistent with available orthogonal results for deletion and copy-number status. At the CYP21A2 locus, TruPath reconstructed the RCCX gene arrangement and assigned the I2G variant and a chimeric allele to different haplotypes. Conclusions: In this small series of archived clinical DNA specimens, TruPath Genome enabled characterization of selected repeat expansions and homologous or pseudogene-associated variants. The findings support further evaluation of proximity-mapped sequencing as a complementary approach for complex genomic regions, particularly when conventional short-read sequencing is limited.

Myotonic dystrophy type 1 (DM1) is caused by a CTG expansion in the 3' untranslated region of the DMPK gene that results in the expression of expanded CUG RNA which trap splicing factors resulting in widespread alternative splicing dysregulation. A panel of these splicing events, or splicing index (SI), is currently used to measure splicing dysregulation in DM1 patients. Despite significant interest in splicing, little is known about its progression over time in DM1. The objective of this study was to evaluate the changes in SI over 3 years in a well-phenotyped cohort of DM1 participants and identify factors contributing to these changes. RNA was extracted from vastus lateralis (VL) biopsies of 20 DM1 participants taken 3 years apart and the SI score was measured using targeted RNA sequencing. Baseline SI scores correlated with percent predicted muscle strength of the ankle dorsiflexors, knee extensors and hand grip. On average, SI scores worsened over 3 years, but the progression varied markedly among participants. Adult/juvenile phenotypes and lower baseline SI scores were among the factors associated with higher progression of SI scores over time. These findings emphasize the importance of integrating clinical phenotype and baseline molecular severity when interpreting longitudinal SI change, particularly in the context of clinical trials.

Pathogenic variants in the myocilin (MYOC) gene are the most common cause of Mendelian open-angle glaucoma. In 2022, the Clinical Genome Resource (ClinGen) Glaucoma Variant Curation Expert Panel (VCEP) published rule specifications for MYOC variant interpretation, including a pilot study of 81 variants. Here, we present the results of curating 271 MYOC variants reported in people with open-angle glaucoma using updated specification rules. Of all the variants, 11 were classified as benign (B), 45 as likely benign (LB), 166 as variants of uncertain significance (VUS), 35 as likely pathogenic (LP), and 14 as pathogenic (P). All LP/P variants were located within the conserved olfactomedin domain encoded by exon 3. The updated variant curation guidelines from the Glaucoma VCEP increased the number of clinically definitive classifications from 28% (74/265) to 39% (105/271), with 95% (41/43) of reclassified variants moving to greater clinical relevance. Functional evidence was lacking for 93% (154/166) of VUS. Additional functional evidence could further enhance classification by halving (85/166) the proportion of those classified as VUS. These findings highlight the role of rule calibration and rigorous functional evidence assessment toward improving variant classification with clinical utility for patients.

Background: Dysfunction of the KCND3-encoded Kv4.3 channel, which carries the transient outward potassium current (Ito), has long been suspected to contribute to the pathophysiology of Brugada Syndrome (BrS). Although recent genome-wide association studies (GWAS) have implicated KCND3, the role of rare variation in this gene remains uncertain and clinically unactionable. Methods: Genome sequencing was performed in 983 BrS probands from Europe and Thailand. Putative gain-of-function (GoF) rare variants in KCND3 identified in patients were functionally characterized by patch-clamp analysis of Kv4.3 current properties using a HEK293 cell heterozygous model. Biobank datasets were utilized to evaluate rare-common variant interaction at the KCND3 locus. Results: Eight BrS cases carried six different rare KCND3 variants which clustered at the Kv4.3 C-terminus, with terminal exon rare variants significantly enriched in cases versus study (p=0.004) and population (p=2.0e-6) controls. Four of the five KCND3 case variants studied demonstrated significant GoF effects versus wild-type and putative benign variants, including increased Ito amplitude (p.Leu586Val, p.Ala634Thr) and slower inactivation (p.Gly600Arg, p.Ala654Thr, p.Ala634Thr). Analysis of p.Gly600Arg carriers in the UK Biobank and All of Us revealed a significant enrichment of the protective KCND3 expression-reducing allele of the GWAS lead single-nucleotide polymorphism rs12090194, suggesting a possible epistatic common-rare variant interaction modifies penetrance in carriers. Conclusions: GoF rare variants at the C-terminus of KCND3 are associated with BrS in a small proportion of cases. Determination of the rare and common variant status in BrS patients and their relatives may provide information about disease risk.

Targeted precision therapies are increasingly used in the treatment of individuals with vascular anomalies (VAs). This increases the need for rapid, accurate and inexpensive genetic diagnosis. Droplet digital polymerase chain reaction (ddPCR) is an alternative to next-generation sequencing (NGS), permitting rapid, highly sensitive interrogation of recurrent pathogenic mosaic variants. We examined the feasibility of ddPCR as a primary diagnostic tool in a large cohort of individuals with VAs. Lesional tissue was collected for ddPCR of up to 46 recurrent pathogenic variants across 16 genes associated with VAs. Specimens were assessed on a subset of assays for each individual based on clinical phenotype. Most individuals who had negative ddPCR results went on to high-depth gene panel or deep exome NGS, or Sanger sequencing. Here we report the phenotypic and molecular findings for 78 newly recruited and tested individuals in addition to the 60 individuals already reported from our cohort. The overall diagnostic yield for our cohort when combined with individuals previously reported was 104/138 (75%). Of 138 individuals tested, recurrent pathogenic variants were detected in 71 (51%) on ddPCR. Variants were most frequently identified in PIK3CA (n=28), TEK (n=18), GNAQ (n=12), or MAP2K1 (n=7). In a further 33 individuals, pathogenic variants were identified on NGS or Sanger sequencing. Our findings indicate that ddPCR is an efficient method achieving a high diagnostic yield in our cohort when used prior to sequencing.

Purpose A significant proportion of cases with rare inherited retinal disease (IRD) remain genetically unresolved following short-read whole genome sequencing (WGS). This study aimed to increase the diagnostic yield of two cohorts in a research setting. One cohort was previously screened by short-read WGS (n=120 cases), while the second cohort (n=28 cases) was not included in any previous short-read WGS studies. In contrast to the larger cohort, these cases only underwent either exome sequencing (ES) or single-molecule molecular inversion probe (smMIPs) sequencing as pre-screening in earlier studies. For all cases, we performed a stepwise, case by case reanalysis in which short-read WGS was generated (n=28 cases) and (re)analyzed short-read WGS data using an optimized approach. For a subset of cases (n=20) that remained unresolved, this was followed by long read WGS. Together, these steps aimed to improve the diagnostic yield in these cohorts. Methods Short-read WGS data of 148 IRD-cases were (re)analyzed using updated allele frequency databases, new variant caller and variant predictor tools, as well as updated and extended gene-panels to identify causal single nucleotide variants (SNVs) and structural variants (SVs). For 20 genetically unresolved cases with sufficient high-quality DNA available, long-read WGS was performed. Results The (re)analysis of short-read WGS data in the 120 previously WGS-screened unresolved cases resulted in a diagnostic yield of 15% (18/120). For the 28 individuals without prior WGS, performing WGS analysis within this study provided an additional yield of 32% (9/28). Furthermore, long-read WGS contributed two genetic diagnoses among the 20 long-read WGS-sequenced probands. Amongst other reasons, new genetic diagnoses were attributed to pathogenic variants in genes newly associated with IRDs and in genes that were not included in the gene panels applied in previous WGS studies. The Mobile Element Locator Tool facilitated the identification of two pathogenic Alu insertions in short-read WGS. Long-read WGS identified two additional pathogenic variants. Conclusion These results support periodic reanalysis of existing short-read WGS data as a key component of IRD genetic diagnostics, with selective long-read WGS representing a valuable additional approach for unresolved cases. Together, these complementary strategies offer a practical framework for improving molecular diagnosis and narrowing the remaining diagnostic gap in IRDs.

Background: Recent work has brought large language models (LLMs) into clinical variant interpretation. In most existing approaches, the model is confined to text-level operations, retrieving, extracting, compacting, and summarizing heterogeneous evidence, while the genuine genetic evaluation is delegated to a separate deterministic engine. Other tools employ LLMs to mine functional evidence for a single ACMG criterion, to rerank a classifier's output, or to answer isolated variant questions. Integrative genetic reasoning, by contrast, demands multi-source evidence weighing for each variant, reconciliation of inheritance with observed zygosity and parental genotypes, joint interpretation of variants in the same gene, and the derivation of a case-level conclusion. Moreover, to protect patient privacy, this must be done on hospital premises. Methods: o build an end-to-end variant interpreter that runs on local infrastructure, we benchmarked 13 small language models (SLMs) on 150 multiple-choice genetics questions and selected Qwen3.5-9B. In IVA (Intelligent Variant Analyser), a four-stage pipeline, the model assigns ACMG/AMP criteria, applies de novo, dominant, recessive, and X-linked inheritance logic to trios, and writes a structured report, while deterministic code handles retrieval, phasing, point summation, citation checking, and formatting. We evaluated IVA on synthetic scenarios, 16 solved trios, and 327 solved singleton cases of visual impairment, on a single GPU. Results: On the benchmark, Qwen3.5-9B reached a mean accuracy of 0.81, compared with 0.89 for a 1.6-trillion-parameter cloud model. In the 16 trios, IVA reported every previously identified pathogenic, likely pathogenic, or uncertain-significance variant of interest (sensitivity 100%; precision 39.5%). In the 327 singleton cases, it recovered 93% of the expert-reported causative variants and classified 82% of them as pathogenic or likely pathogenic (precision 43%). Failures fell into four recurrent classes: omission, unverified mechanism, mismatch between reasoning and label, and arithmetic errors. Conclusions: A small language model running on local hardware can serve as the reasoning layer for case-level and trio-based variant interpretation when it is surrounded by deterministic safeguards. Because its outputs are not deterministic, analyst review and prospective multicenter validation remain necessary.

Centronuclear myopathy (CNM) is a rare form of inherited diseases often caused by single base mutations. Modelling CNM is challenging due to the diversity of CNM mutations and genetic background of individual patient. To address this, we used base editing to introduce CNM mutations into an induced pluripotent stem cell (iPSC) line from a healthy donor, generating a collection of iPSC lines (iPS-CNM) carrying distinct CNM mutations. We found that the efficiency of base editing depended critically on selecting base editor (BE) variants and the target sequences. Optimization using different BE variant and sgRNA pair was required for each target site. Moreover, whole genome sequencing (WGS) was performed to confirm on-targets and detect off-targets in order to select iPSC clones for the collection. Our findings highlight the feasibility of base editing for generating an iPSC collection from one parental iPSC line combined with thorough evaluation of rare off-targets using WGS.

Background: Hereditary spastic paraplegia (HSP) comprises a heterogeneous group of neurodegenerative disorders characterized by progressive lower-extremity weakness and spasticity, with or without additional neurological symptoms. Objectives: To characterize the molecular genetic spectrum, inheritance patterns, molecular finding rate, and clinical characteristics of HSP within a multicenter Canadian cohort. Methods: We investigated 671 patients from 423 families who met clinical criteria for HSP identified through the Canadian HSP (CanHSP) registry from 7 medical centers across Canada. Patients were evaluated using known HSP and HSP mimic gene panels, whereas unsolved and more recent patients underwent whole-exome sequencing (WES). 435 individuals underwent gene-panel testing and 485 underwent WES. Copy-number variants (CNVs) were assessed using ExomeDepth, and candidate CNVs using multiplex ligation-dependent probe amplification (MLPA). Results: Molecular findings were identified in 426 of 671 individuals (63.5%): 351 (52.3%) had HSP-associated findings, 70 (10.4%) HSP-mimicking disorders, and five (0.7%) other molecular findings; 245 (36.5%) remained unresolved. Candidate variants were found in 80 relevant genes, 36 of which were HSP-associated genes. The most frequently implicated genes were SPAST (22.5%), SPG7 (5.7%), ATL1 (5.2%), SPG11 (3.6%), and SPG8 (3.3%). Phenotypically, SPAST, ATL1, and SPG8 w[zg1.1]ere predominantly associated with pure HSP, whereas SPG7 and SPG11 more often showed complicated phenotypes, with prominent cerebellar involvement in SPG7 and cognitive and developmental features in SPG11. Conclusion: Patients in CanHSP present with genetic and phenotypic heterogeneity. Reanalysis and complementary approaches, including whole-genome, RNA, and long-read sequencing, may increase diagnostic yield.

Using a genetic screening approach based on an inducible gene-activating system and cell sorting, Down syndrome critical region 3 (DSCR3) was isolated as a gene whose overexpression increased cell size. Fibroblasts derived from individuals with Down syndrome (DS) exhibit elevated DSCR3 expression at both the mRNA and protein levels, correlating with increased cell volume compared to fibroblasts from healthy donors. Despite a slower proliferation rate, DS fibroblasts demonstrate higher basal and maximal mitochondrial respiration, suggesting enhanced metabolic activity associated with increased cell size. siRNA-mediated knockdown of DSCR3 reduces cell size in both DS and normal fibroblasts, indicating its general role in cell size regulation. As DSCR3 is a component of the retriever complex involved in endosomal cargo recycling, these findings position membrane protein trafficking as a novel module for cell size control.